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Abstract 


In this paper we study the robustness with respect to stability of the 
closed-loop system with collocated rate sensors using LQG (mean square rate) 
optimized compensators. Our main result is that the transmission zeros of the 
compensator are precisely the structure modes when the actuator/sensor 
locations are “pinned” and/or “clamped”: i.e., motion in the direction sensed is 
not allowed. We have stability even under parameter mismatch, except in the 
unlikely situation where such a mode frequency of the assumed system 
coincides with an undamped mode frequency of the real system and the 
corresponding mode shape is an eigenvector of the compensator transfer 
function matrix at that frequency. For a truncated modal model — such as that 
of the NASA LaRC Phase Zero Evolutionary model — the transmission zeros 
of the corresponding compensator transfer function can be interpreted as the 
structure modes when motion in the directions sensed is prohibited. 
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1. Introduction 


The robustness of control laws to parameter uncertainty is of particular importance 
to Space applications because testing large structures under the micro-g conditions in Space 
is not possible on the ground. This paper explores the robustness issue for LQG optimized 
compensators using the explicit form discovered by the author for their time and/or 
frequency domain representation, whatever the structure model used — whether it is FEM, 
Truncated Modal, or Continuum. 

The basic properties affecting robustness of the LQG optimized compensator are 
developed in Section 2 and how they relate to robustness is examined in Section 3. 
Section 4 deals with continuum models where in particular we show that the transmission 
zeros are the modes of the structure when the actuator/sensor locations are “pinned” and/or 
“clamped,” i.e., motion in the directions sensed is restricted, generalizing the usual 
notions for simple beams. 

When a truncated modal model is available, the compensator can be expressed 
explicitly also in terms of the given modes and mode shape vectors. The transmission 
zeros of these approximate compensators are studied in Section 5 and in particular some 
numerical results are presented for the NASA LaRC Phase Zero Evolutionary Model. 
Conclusions are in Section 6. 


2. The LQG Optimized Compensator 

To state the LQG problem, we begin with the canonical time-domain dynamics of a 
flexible structure with collocated rate sensors which, whether it is a Finite Element Model 
or Truncated Modal model (and hence finite dimensional) or a Continuum Model (and 
hence infinite-dimensional), can be expressed in the form: 


Mx{t) + Ax(t) + Bu{t) + BN a (t) = 0 
v(r) = B*x(t) + N r (t) 


(2.1) 


where in the case of FEM, 

M is the mass matrix (nonsingular, nonnegative definite) 

A is the stiffness matrix (nonsingular, nonnegative definite) 

B is the control matrix (rectangular matrix) 

w(-) is the control vector («xl, assuming n actuators) 

*(•) is the “displacement” vector 

N a (-) is the actuator noise assumed white Gaussian with spectral 

density d Q l , / being the nxn Identity matrix 

v(*) is the sensor output 

B* represents the transpose of B 

N r (-) is the sensor noise assumed white Gaussian with spectral density d r I . 



For the Continuum Model such a representation continues to hold, however compli- 
cated the structure, with *(•) now allowed to range in a Hilbert space X, with A, M, B 
being linear operators: 

M bounded linear, self-adjoint, nonnegative definite with M 1 bounded; 

A closed linear, self-adjoint, nonnegative definite with compact resolvent, 

the resolvent set including zero 

B maps E" Euclidean n-space into X, and 

B* represents the adjoint of B. 

See [1, 2], 

The LQG problem we shall consider is that of finding the control «(•) (or equiva- 
lently the optimal compensator) that minimizes the mean square time average of the rate: 

lim Uf \\B*x(t)\\ 2 dt + h J ||w(f)|| 2 dt\ (2.2) 

T-- l' o o J 

where X. > 0. 

It is shown in [1, 3] that under the “controllability” assumption that 

B*§ k * 0 

for any k, where <j> 4 are the modes orthonormalized with respect to the mass matrix: 

i4<j)* = co^Aftt)* ; [M$ k , <)>*] = 1 , (2.3) 

the optimal compensator transfer function (nxn matrix function) can be expressed in the 
explicit analytical form: 

\\t(p) = gpB*(p 2 M + A + ypBB*) 'fi , Re. p > 0 (2.4) 


where 


g = 


<djd r 


Y . + i . 


Moreover, the corresponding mean square control power is given by 

T d„ 


lim \\ \m\\ 2 dt = Tr. (B*MB )~ 1 


0 

and the corresponding mean square displacement is: 

T JKd 

||B**(0ll 2 dt = <djd r + 

o 


(2.5) 


(2.6) 


. T ( VX d 

lim =J ||B*x(0ll 2 dt = ^d a /d r + — j Tr. BM B . (2.7) 

T -4 00 n ' 


See [4] for the corresponding time-domain version of (2.4). From (2.3) we can 
deduce readily that 



(i) As X-*0, g -i °° and y -> °° we note that ^ -> 'Sdjd r and hence 


\\f(p) -> Ndjd r 

where the right side is recognized as is the optimal “static” or “direct connection” or 
“PID” controller. Note that as X -> 0, the control power given by (2.6) becomes infinite, 
as we expect. 

(ii) \\r(p) is “positive real” — that is to say: 

\|/(p) holomorphic in R e.p > 0 

vj/(p) + \|/(/?)* nonsingular, and positive definite, for Re. p > 0 

where * denotes conjugate transpose. We shall prove this directly, here, even though it 
may be deduced from the results in [1], In fact 

( p 2 M + A + ypBB*)x = 0 


implies that 


p 2 [Mx, x] + [Ax, x] + yp\\B*x\\ 2 = 0 


and normalizing so that 


[Mx, x] = 1 


we obtain 

p 2 + yp\\B*x\? + [Ax, x] = 0 . 

Because of our assumption that B*§ k is not zero for any k, we see that 


115**11 > 0 


and hence 

Re. p < 0 . 

This is enough to imply that in the finite-dimensional case: 

(p 2 M + A + y pBB*)~ l 

is holomorphic in Re,/? > 0. In the infinite dimensional case the fact that A has a 
compact resolvent implies that so does 

p 2 M + A + y pBB* 

and hence it follows that 

B*(p 2 M + A + y pBB*)~ l B 

is holomorphic in Re.p > 0, and in fact is an function. 

Next let us calculate 

V(P) + V(P)* > Re.p >0 . 


448 



We have 


¥(P) ± VO)* _ B*(p 2 M + A + ypBB*) 1 

• [|/?| 2 (/? + p)M + (p + p)A + 2y\p\ 2 BB*](p 2 M + A + y pBB*) 1 B 


which is >0, since 

\P\ 2 (P+P)M + ( P+P)A + 2y\ p\ 2 BB* 


is, for Re. p > 0 . 

In particular for 


p = ic o , 


-oo < (0 < OO 


we have 

±. ¥( t(a ) _ y(0 2 \)/((C0)V|f(/C0)* . (2.8) 

This leads to an important result which we state as: 

Lemma 2.1 

Suppose for some (0, < © < «, 

Re. [\|/(/co)v, v] = 0 for some v . (2.9) 

Then 

\|/(i (0)v = 0 (2.10) 

and if co * 0, 

B*(-co 2 Af + /4) -1 Rv = 0. (2.11) 


Proof. 


which by (2.8) is 


Re. [V|/(i"co)v, v] 


[( 


M/Q'cd) + \|/(/co)* 


*vll 2 


)v. v] 


and hence (2.9) is equivalent to: 


= Yar||\j/(iCD)*v 


co||v|/(j'©)*v|| = 0 . 


If (0 = 0, \|/(fco) = 0 and hence y(/(0)v = 0. 

If (0 is not zero, 

y(/co)*v = 0 


or 


Let 


or 


6 *(-co 2 M + A - iycpBB*) 1 Bv = 0 
( — co 2 A/ + A - iyaiBB*) 1 Bv = x 
Bv = -vfMx + Ax - iyozBB*x . 
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Since 


B*x = 0 , 


we have 


B\ = ( -G) 2 A/ + A)x 


and co is not an eigenvalue of A, since B*x = 0. Hence 

B*(-co 2 A/ + A )~ 1 Bv - 0. 


Corollary. 

Except for co = 0, the zeros of 


Det |v|/(/co) + \|/(jco)*| 

in -oo < co < « are the same as those of 

Det |B*(-co 2 A/ + .<4) — 1 J5| = 0 
and in particular independent of g and y, for g + y < °°. 


Proof. 


implies that 


and hence by the lemma, 


and since co * 0, 


Det |y(<co) + vj/(ico)*| = 0 


[(vj/(i'co) + \|/(ico)*)v, v] = 0 


\|/(ko)*v = 0 


B*(-co 2 M + A) Bv = 0 


Conversely if 


Det lB*(-co 2 Af + A) = 0 . 


Det |fi*(-co 2 M + A) 1 S| = 0 , 


so that 


(B*(-C0 2 M + A) X B)v = 0 


(-0 ) 2 M + i4) B\ = x 


Bv = -co 2 Mjc + Ax = -CO 2 Mx + Ay ± y(iu>)BB*x 
for arbitrary value of y, since B*x = 0. Hence 

B*(-(0 2 M + A ±y(/co)BB*) _1 Bv = 0 
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or 


(V(/CO) + Y(i'CO)*)v = 0 
or 

Det |\|/(i'co) + y(('co)*| = 0 . 


3. Robustness 


The robustness of concern is that with respect to parameter uncertainty, in particular 
in the mode frequencies CD*. Thus we want to be able to assert that the closed-loop 
system is stable even if the parameters chosen for the compensator transfer function \|/( ) 
are incorrect. Now the closed-loop transfer function corresponding to the compensator 
transfer function y(-) is given by 

(p 2 M + A + ypB\y(p)B*y l . (3.1) 

Let p be a pole of (3.1), so that for some x * 0, 

p 2 Mx + Ax + ypBy(p)B*x = 0 . (3.2) 


Then 


p 2 [Mx, *] + [Ax, x] + yp[B\\f(p)B*x, jc] = 0 


where [Mx, x] > 0 . 
If p = 0, then 


Ax = 0 


which is not possible since zero is not an eigenvalue of A — this property of the system 
is assumed to be known with certainty. If p * 0, we may divide through by p to get 


Let 


p[Mx, x] + \- Ax ’ 

P 


+ y [B\\i(p)B*x, x] = 0 . 


p = a + /co , 


a > 0 . 


Then 


a [Mx, x] + ^^’^2 + Y[ B (V(P) + y(p)*)B*x, x] = 0 . (3.3) 

Suppose a > 0. Then by the positive real property, the third term in (3.3) is positive and 
the first two terms are of course positive, and hence the sum cannot be zero and hence a 
cannot be positive. Consider next the case a = 0. This yields 

(ia)[Mx, x] - 1 + y[B\\t(i(ti)B*x, *1 = 0. 

For this to hold, it is necessary that 

Re. x\ = 0 

where B*x cannot be zero. Hence, by the lemma 

= 0 , B*x * 0 . 
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B*x * 0 


(3.4) 


But this in turn by (3.2) would imply that 

-m 2 Mx + Ax = 0 ; 

\y(i(j))B*x = 0 

In other words an undamped mode-frequency of the system must coincide with a zero of 
\|/(i(o) corresponding to the same value of B*x (“mode shape at the sensor location”). 
Thus we have robustness with respect to stability so long as this is insured against. Be- 
cause of the coincidence requirement on the mode shape in addition to the frequency, this 
is highly unlikely if the controller dimension (number of actuators) is higher than one. 


4. Continuum Models 

For the case of the continuum model, whether explicit or conceptual, we can relate 
the zeros of the compensator transfer function to the model in a simple way — viz., we 
can show that they are the “pinned” and/or “clamped” mode frequencies of the structure 
or a slight generalization thereof. 

We begin with the general case of a multi-beam model as the NASA LaRC Phase 
Zero Evolutionary Model [5], Here the state variable jc(*) has the form (see [2]): 

/ 

x - 

b 

where /(•) represents the displacement (6x1) vector and b the corresponding “boundary” 
values at the nodes, and thus a finite dimensional vector. Also 

M of 

Mx = 

M b b 
Aof 

Ax = x e 3) (A) 

A bf 


Bu = _ (where B*B U is nonsingular) 

t> u u 

where the dimension of U (the control vector) can be smaller than that of b, and B u 
maps U into the finite-dimensional space spanned by b. In this case 

Det |B*(-t0 2 A/ + A)~ l B\ = 0 (4.1) 

and equivalently, for some v : 

(-(0 2 Af + A)x = B\ ; B*x = 0 ; 

which under the notation: 

f 

x = 

b 
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becomes: 


-m 2 M 0 / + A 0 f = 0 


(4.2) 

(4.3) 

(4.4) 


B*b = 0 

-CO 2 M b b + A*/ = 

The condition (4.4) is superfluous since v is not specified and all we require is that v 
and hence B u \ be nonzero. Thus the transmission zeros are the eigenvalues of the 
differential operator A 0 with the “boundary condition” specified by 


-<o 2 Mo/ + Aof = 0 ; B*b = 0 . (4.5) 

These are recognized as structure modes when the control/sensor locations are “pinned” 
and/or “clamped” — motion in the directions sensed is not allowed. The structure modes 
when all nodes are clamped are of course given by: 


-co 2 M 0 / + A 0 f = 0 , /* 0 


b = 0 . 


(4.6) 


If all nodes are control nodes so that B u is the identity, these are also transmission zeros 
but not in general because of the additional condition (4.4). 

As shown in [3], (4.1) can be further reduced to: 

|B*(-0>X + 1Xto))~ l B H \ = 0 (4.7) 

where the “clamped” mode frequencies given by (4.6) are the “poles” of the matrix T(ia>) . 

A textbook example of (4.7) is provided by the torsion of a one-dimensional beam 
with one end clamped and the other end the control node. Here T(/(o) is given by (see 
[3] for details): 

T(ic o) = VpG / y to cot (2£Vp/G )oa 
and hence the compensator transmission zeros are given by 

sin (2£VpA3)co = 0 


or 


whereas the structure undamped modes are given by 

mco - (Vp G)/^ cot (2£Vp/G)co . 

Note that asymptotically these frequencies merge — a phenomenon which can be proved 
to hold generally. In this (one-dimensional) case the zeros and poles of 

(-0 ) 2 m + T(/co)) 

alternate (an instance of Foster's Theorem familiar in classical circuit analysis [6]) and 
hence also the compensator transmission zeros and undamped mode frequencies — but 
this is no longer true in general in the multidimensional case. 




5. Truncated Modal Model 


Truncated Modal models provide both the (undamped) mode frequencies 00 * and the 
corresponding mode shapes (column vectors of dimension equal to the control dimension) 
B*$ k at the sensor locations up to a maximum frequency allegedly adequate for faithful 
representation of the structure. Assume thus that (0 t , B*§ k are given for k = 1, N. 
Then we may consider this the “truth model” and the corresponding optimal compensator 
transfer function takes the form (see [4]): 

8P B $(p 2 l + D N + "lP B N B ftY lB N ( 51 ) 

where 

D n = Diag. ((01, (Djv) 

Bjf - I B*< (>, B *( |> 2 • • • B*< \ nxN 


B n = 


(**<!> t)* 

( b *$n )* 


Nxn 


where n is the control dimension (equivalently, the number of actuators). It is assumed 
that N is large enough so that the nxn matrix: 


B8 b n 

is nonsingular. Then omitting co = 0, the transmission zeros of the transfer function (5.1) 
are given by 

Det |B$(-(0 2 / + D n ) l B N | = 0. (5.2) 


For large enough N we should expect these frequencies to closely approximate the 
structure mode frequencies when motion is restricted along the directions sensed at the 
actuator/sensor locations. In particular if the theoretical values of the latter are known, we 
have a means of checking the faithfulness of the truncated model. 

For the NASA LaRC Phase Zero Evolutionary Model [5], the truncated modal model 
has 86 modes. For the corresponding mode shapes as determined by the LaRC team, the 
frequencies for which (5.2) hold are shown in Figures 1 through 6, where the minimum 
absolute value of the eigenvalues of 

Bfi (—CO 2 / + D n )~'b n 


is plotted as a function of omega for N = 86. Note that all the eigenvalues are positive 
for a) < (0i and negative for <0 > Figure 1 shows the entire range from 0-300 
radians/second. Figures 2-6 show more detail of the behavior over narrower ranges. 

The dependence of zeros on the depth of the modal approximation is illustrated in 
Figures 7, 8, 9, la, 2a, 3a, 4a, 6a and 10a for N = 8, 16, 30, corresponding to 
( 0 8 = 10.921, (0 16 = 25.225, 0)30 = 53.132, respectively. Note that for the 8-mode model 
there are no zeros at all, while the 16-mode model shows three zeros (in the range 0-25 
rad/sec). The 30-mode model shows excellent agreement with the 86-mode model for (0 
up to (030 = 53.132, comparing Figures 1 and la, 2 and 2a, 3 and 3a, 4 and 4a, 6 and 6a, 
and finally Figures 10 and 10a show the expected divergence for (0 > ( 030 - For illustrative 
purposes we list the first few zeros for the full 86-mode model in rad/sec: 
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(5.4913) 

6.575 

8.75 

(9.2580) 

(14.46) 

14.7 

(15.26) 

where the numbers in parentheses are the nearest undamped mode frequencies. 


6. Conclusions 

It is shown that the optimal compensator transfer function for LQG rate minimization 
for flexible structures with collocated rate sensors has transmission zeros at frequencies to 
the modes of the corresponding continuum structure when the control-sensor locations are 
“pinned” and/or “clamped” (motion in the directions sensed is curbed). In particular the 
compensator is robust with respect to stability so long as any such mode of the assumed 
system does not coincide with an undamped mode frequency of the real system and the 
corresponding mode shape at the sensor locations is an eigenvector of the compensator 
transfer function matrix at that frequency. For Continuum Models the transmission zeros 
are shown to be the poles of a matrix function related to the undamped modes. Calcu- 
lations of the zeros are given for the truncated modal models of the NASA LaRC Phase 
Zero Evolutionary Model illustrating the dependence on the number of modes used. 
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Figure 5 

Minimum Absolute Eigenvalue: 86-mode Truth Model 
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Figure 6a 

Minimum Absolute Eigenvalue: 30-mode Truth Model 
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Figure 7 

Minimum Absolute Eigenvalue: 8-mode Truth Model 
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